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PACS 87. 16. Ka - Filaments, microtubules, their networks, and supramolecular assemblies 
PACS 87. 10. -e - General theory and mathematical aspects 

Abstract - We study analytically and numerically a generic continuum model of an isotropic 
active solid with internal stresses generated by non-equilibrium 'active' mechano-chemical reac- 
tions. Our analysis shows that the gel can be tuned through three classes of dynamical states by 
increasing motor activity: a constant unstrained state of homogeneous density, a state where the 
local density exhibits sustained oscillations, and a steady-state which is spontaneously contracted, 
with a uniform mean density. 
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Introduction. — The mechanical properties of living 
cells are largely controlled by a variety of filament-motor 
networks optimized for diverse physiological processes [l]. 
In the presence of ATP, such networks are capable of gen- 
erating controlled contractile forces and spontaneous os- 
cillations. These phenomena arise from the presence of 
groups of motor proteins, such as myosin II, that convert 
the chemical energy from ATP hydrolysis into mechanical 
work via a cyclic process of attachment and detachment to 
associated polar protein filaments, e.g. F-actin [2]. There 
is great variation in the organization of actin, myosins, 
and other cross-linking proteins in the actomyosin struc- 
tures found in cells. Myofibrils in striated muscle cells are 
examples of highly organized structures [Tj , composed of 
repeated subunits of actin and myosin, known as sarcom- 
eres, arranged in series. Each sarcomere consists of actin 
filament of alternating polarity bound at their pointed end 
by large clusters of myosins, known as myosin "thick fila- 
ments". The periodic structure of the myofibril allows it 
to generate forces on large length scales due to the collec- 
tive dynamics of individual units of microscopic size, giv- 
ing rise to muscle oscillation and contraction. More diffi- 
cult is to understand the origin of spontaneous oscillations 
and contractility in cytoskeletal filament-motor assemblies 
that lack such a highly organized structure [3] . 

Oscillations are for instance observed in many organ- 
isms during repositioning of the mitotic spindle from the 
cell center towards the cell pole when unequal cell divi- 



sion occurs [3J. In vitro examples of such phenomena are 
sustained cilia-like beatings in self-assembled bundles of 
microtubules and dyneins [I] and spontaneous contractil- 
ity in isotropic reconstituted actomyosin networks with 
additional F-actin crosslinking proteins, such as filamin or 
a-actinin 151161 . 



Theoretical models have shown that spontaneous os- 
cillations can arise from the collective action of groups 
of molecular motors coupled to a single elastic ele- 
ment [7] . In these models the load dependence of the bind- 
ing/unbinding kinetics of motor proteins breaks detailed 
balance and provides the crucial nonlinearity that tunes 
the system through a Hopf bifurcation to a spontaneously 
oscillating state. This simple theoretical concept has been 
adapted to describe the beating of cilia and flagella |8j, 
mitotic spindles during asymmetric cell division [9], and 
spontaneous waves in muscle sarcomeres 10 -131. 



In a parallel development, generic continuum theories 
of active polar 'gels' have been constructed by suitable 
modification of the hydrodynamic equations of equilib- 
rium liquid crystals to incorporate the effect of activity. 
These continuum models are capable of capturing some 
of the large-scale consequences of the internal contractile 
stresses induced by active myosin crosslinkers, including 
the presence of propagating actin waves in cells adher- 
ing to a substrate 14 and the retrograde flow in the 
lamellipodium of crawling cells 15 . In these studies the 



acto-myosin network is modeled as a Maxwell viscoelastic 
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fluid, with short-time elasticity and liquid response at long 
times |16| . The active viscoelastic liquid cannot, however, 
support sustained oscillations that require low frequency, 
long wavelength elastic restoring forces. 

In this letter we consider a nonlinear version of the 
generic continuum model of an isotropic active solid dis- 
cussed recently by two of us 17 . We go beyond the linear 



stability analysis discussed in 17 and show that the pres- 
ence of nonlinearities leads to both stable contracted and 
oscillatory states in different regions of parameters. The 
acto-myosin network is modeled as an elastic continuum, 
as appropriate for a cross-linked polymer gel embedded in 
a permeating viscous fluid, with elastic response at long 
times and liquid-like dissipation at short times. The active 
solid model describes the various phases of acto-myosin 
systems as a function of motor activity, including spon- 
taneous contractility and oscillations. It provides a uni- 
fied description of both phenomena and a minimal model 
relevant to many biological systems with motor-filament 
assemblies that behave as solids at low frequencies. When 
the active solid is isotropic, as assumed in the present 
work, the coupling to a non-hydrodynamic mode provided 
by the binding and unbinding kinetics of motor proteins 
is essential to generate spontaneous sustained oscillations. 

The main results are summarized in Fig. [TJ where we 
sketch the steady states of the system as we tune the ac- 
tivity, defined as the difference Ap between the chemical 
potential of ATP and its hydrolysis products, and the com- 
pressional modulus B of the passive gel. For a fixed value 
of B we find a regime where the active gel supports sus- 
tained oscillating states and a contracted steady state as 
Ap is increased. For fixed Ap spontaneous contractility 
is only observed below a critical network stiffness. This is 
consistent with experiments on isotropic acto-myosin net- 
works with additional cross-linking a-actinin where spon- 
taneous contractility was seen only in an intermediate 
range of a-actinin concentration [5J. Our model does not, 
however, yield a lower bound on B below which the con- 
tracted state does not exist. This may be because, in con- 
trast to the experiments where a minimum concentration 
of a-actinin is required to provide integrity to the net- 
work, our system is always by definition an clastic solid, 
even at the lowest values of B. The phase diagram result- 
ing from our model also resembles the state space diagram 
of a muscle sarcomere obtained experimentally 
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Interestingly, we find that a simplified dynamical system 
obtained by a one-mode approximation to our continuum 
theory corresponds to the half- sarcomere model proposed 
recently by Giinther and Kruse 10 for a particular set of 
parameters. Our analysis shows, however, that the phase 
behavior described above is generic and can be expected in 
a wide variety of active elastic systems, as it relies solely on 
symmetry arguments. It provides a unified description of 
both spontaneous and oscillatory states and predicts their 
region of stability as a function of the elastic properties of 
the network and motor activity. 
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Fig. 1: (Color online) "Phases" of the nonlinear active elastic 
gel obtained by varying the compressional modulus B and the 
activity Ap. 



Model. - In Ref. 17 



two of us formulated, on the 
basis of symmetry arguments, a phenomenological hydro- 
dynamic model of a cross-linked gel (a network of actin fil- 
aments crosslinked by filamins or other "passive" linkers) 
under the influence of active forces exerted by clusters of 
crosslinking motor proteins (e.g., myosin II minifilaments). 
Only linear terms were retained in the continuum equa- 
tions in [l7| where the linear modes of the system and 
their stability were analyzed. Here we consider a non- 
linear continuum model and show how nonlinearities can 
stabilize oscillatory and contracted states. The active gel 
consists of a polymer network in a permeating fluid |17| . 
On length scales large compared to the network mesh size, 
the gel can be described as a continuum elastic medium, 
viscously coupled to a Newtonian fluid. The model follows 
closely that formulated for a passive gel 18 . We focus on 
compressional deformations and consider for simplicity a 
one-dimensional model. We assume the permeating fluid 
to be incompressible and consider the case of small volume 
fraction of the gel. In this limit the permeating fluid sim- 
ply provides a frictional drag to the polymer network. The 
hydrodynamic description is then obtained in terms of two 
conserved fields: the density p(x, t) of the gel and the one- 
dimensional displacement field, u(x,t). The momentum 
density of the gel is not conserved due to drag exerted by 
the permeating fluid. In addition, density variations are 
determined by the local strain, with Sp = p — po = —pod x u 
and po the equilibrium mean density. The elastic free 
energy density of the gel can be expanded in the strain 
s = d x u about the state s — 0. To describe the possibility 
of swelling and collapse of a gel (even a passive one) one 
needs to keep terms up to fourth order in s in the elastic 
free energy density (dropping a linear term that can be 
eliminated by redefining the ground state) [19] 
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where B is the longitudinal compressional modulus of the 
gel and a, f3 > are phenomenological parameters cap- 



p-2 



Generic phases of cross-linked active gels: Relaxation, Oscillation and Contractility 



turing the effects of many-body interactions and nonlinear 
elasticity of the components |20] , 

It is straightforward to show that within mean-field the- 
ory the free energy given in Eq. ([I]) yields a line of first 
order phase transitions at B = B c = 2a 2 / (9/3) between 
an unstrained state with s = for B > B c and a strained 
state with finite s for B < B c . The stable strained state 
is one of higher density (s < 0) for a > 0, corresponding 
to a contracted state, and one of lower density (s > 0) for 
a < 0, corresponding to an expanded state. The transition 
line terminates at a critical point at a = 0. 

Activity is induced by the presence of a concentration 
c(x, t) of bound active cross-linkers that undergo a cyclic 
binding/unbinding transformation fueled by ATP, exert- 
ing forces on the gel. The dynamics of the active gel on 
time scales larger than the Kelvin- Voigt viscoelastic re- 
laxation time is described by coupled equations for u(x, t) 
and c(x,t), given by 

Td t u = d x a e - d x p a , (2) 
d t c = -d x (cd t u) - k(s)(c - c ) (3) 

where a e — — Bs + as 2 + /3s 3 is the elastic stress, T is 
a friction constant describing the coupling of the polymer 
network to the permeating fluid, and p a (p, c) is the active 
contribution to the pressure, describing the isotropic part 
of the active stresses induced by myosins. Equation ^ al- 
lows for convection of bound motors by the gel at speed dtu 
and incorporates the binding/unbinding dynamics, with 
k(s) the strain-dependent motor unbinding rate and Co the 
equilibrium concentration of bound motors. Since highly 
non-processive motor proteins such as myosins are on av- 
erage largely unbound, we neglect the dynamics of free 
motors that provide an infinite motor reservoir. 

The active pressure is taken to be linear in the rate Ap, 
of ATP consumption, with p a — A/z ((p, c). This is a rea- 
sonable approximation for weakly active systems, where 
the number of active elements make up a small fraction 
of the total mass of the gel, as is the case in most experi- 
ments. []] We then expand 



c(p, c) ^Co - &8p~ - C2</> - C3 (spy 

+ C 5 2 + C 6 (<5,5) 3 ... 



(4) 



with 5p = Sp/po and <j> — (c—cq)/co the fluctuations in the 
gel and bound motor concentrations, respectively, and all 
parameters Q defined positive. The positive sign of £i and 
£2 corresponds to a contractile acto-myosin system and de- 
scribes the reduction in the longitudinal stiffness of the gel 
from contractile forces exerted by motor proteins. The pa- 
rameters (3, (4 and Cs describe excluded volume effects. A 
positive £ 3 favors contracted over expanded states. A pos- 
itive ^6 guarantees the stability of the network in regions 
of negative effective compressional modulus. 



x We note that the effects of nonlinearities in A/x can taken ac- 
count of by expanding in 5fi = A/j, — Afig about a stationary state 
with finite AfiQ. 



There are several sources of nonlinearities in Eqs. (2p): 
the nonlinear strain dependence of the elastic free energy, 
the nonlinear terms in the active pressure, and the depen- 
dence of the motor unbinding rate k on the load force on 
bound motors, which in turn is proportional to the strain 
s of the elastic gel. We assume an exponential dependence 
of the form k(s) — k$e e8 21 , with ko the unbinding rate 



of unloaded motors and e a dimensionless parameter de- 
termined by microscopic properties of the motor-filament 
interaction. In the following we expand the unbinding 
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rate for small strain as k(s) ~ fco 
Keeping higher order terms in e does not change qualita- 
tively the behavior described below. Finally, the first term 
on the rhs of Eq. ([3| contains a convective nonlinearity 
that does not affect the key features of dynamics observed 
and will be neglected in most of the following. 

The equations for the active gel can then be written as 

Td t u = 8 x a a e + Apd x [C 2 + Uscj) ~ ( 5 4> 2 ] , (5a) 

J2 



d t (j) = -d x [(1 + <p) d t u] - k 
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<t> , (5b) 



where cr° = B a s + a a s 2 + /3 Q s 3 , with renormalized elastic 
constants, B a = B — CiAMi a a = a + (3AP and /3 a = 
+ CeA M . 

We render our equations dimensionless by letting u — > 
u/L and t — > tk^ 1 . We then define dimensionless param- 

rfc L 3 



eters as B = ^2 



A/5 



- 8 - g 

rfeo-L 2 ' p Tk L 2 ■ 

and Q — QL. In the following we drop the tilde to sim- 
plify notation and all quantities should be understood as 
dimensionless. 

Linear stability analysis. — There are three steady 
state solutions of Eqs. ( 5a|5b I: an unstrained state, with 
(u s ,(f) s ) — (0,0), and two strained states, with (u s ,<p s ) = 

(s±x, 0) and s± = (^-a a ± ^a 2 a ~ 4B a /3 a ) /2/3 a , provided 

B a < a 2 /4/3 a . For concreteness we choose a a ,/3 a > 0. 
Then if B a > 0, s± < and 5p± > 0, so that both strained 
solutions correspond to contracted states. If B a < 0, then 
s_ < and s + > 0, corresponding to contracted and 
expanded states, respectively. 

In this section we examine the linear stability of each of 
these states. Letting s = sq + 5s, where sq = (0, s±) rep- 
resents the constant value of strain in the steady state and 
Ss = d x du, the linearized equations for the displacement 
and motor concentration fluctuations are given by 



where 



d t Su = Bd 2 Ju + Zd x (j) , 
dt<f> = -d x d t 5u - n4> , 



B = B a + 2a a s + ZPasl 
Z = Ap(( 2 + C 4 s ) , 

K = 1 + g.S + ESq/2 ■ 



(6b) 



(7a) 
(7b) 
(7c) 
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Fig. 2: (Color online) Phase diagram obtained by analyzing 
the linear stability of the modes. The left frame is for £2 = 
£4 = 0, corresponding to <f> — 0. The right frame is obtained 
including motor fluctuations, with £2 = C an d C4 — vC- The 
other parameter values are a — 0.1, /3 — 0.1, £1 = £ = 1, 
C3 = rjC, and ( 6 = V 2 (, with V = 0.1. 



Looking for solution of the form 8u, <f> ~ e zt+iqx we g n j 
that the dynamics of fluctuations is controlled by two 
eigenvalues, with dispersion relations 



b(q) 



± 



1 



6(g)] 2 - 4nBq 2 



(8) 



where b(q) = n + (B—Z)q 2 and z u (q) and z < p(q) correspond 
to the root with the plus and minus signs, respectively. 

If motor fluctuations are neglected (0 = 0), correspond- 
ing to letting Z = in Eq. ^8j, one finds z u = —Bq 2 and 
the linear stability of the steady states is entirely deter- 
mined by the sign of B, with B > 0, corresponding to a 
linearly stable state. When <j> — the problem is equiv- 
alent to an equilibrium problem, with q 2 B the curvature 
of a free energy of the form given in Eq. ([IJ, albeit with 
elastic constants renormalized by activity. In this case the 
unstrained state s = is the global minimum of the free 
energy for B a > 2a 2 /(9/3 a ), while the contracted state 
s = s_ is the global minimum for B a < 2a 2 /(9f3 a ). The 
line B a = 2a 2 /(9(3 a ) defines a line of first order phase 
transitions in the (B,Afi) (see Fig. [2j left frame). If, on 
the other hand, we consider the problem dynamically, the 
condition of stability of linear fluctuations given by B > 
is a necessary, but not sufficient one to identify the stable 
steady states of the system, as multiple fixed points with 
different basins of attraction coexist in the same region 
of parameters. In particular, both the unstrained s = 
state and the contracted s = s_ state are linearly stable 
for B a < B < a 2 /(4/3 a ), while both the contracted and 
expanded states (s = s T ) are stable for B a < 0. If we 
assume that among these linearly stable states the system 
selects dynamically the steady state with the fastest de- 
cay rate q 2 B, then we recover the linear phase diagram of 
Fig. [2] obtained form the equilibrium analysis. 

When the coupling to motor concentration fluctua- 
tions is included, we find a region of parameters where 
Re [ z u,<t>] > for all three homogeneous solutions s = 



(0, s±). This is the white region in Fig. [2] (right frame), 
where the dynamical system is linearly unstable. The in- 
stability occurs in a region where the modes are complex, 
describing oscillatory states, and Z — njq 2 > B > 0. The 
linear stability diagram for the system is shown in Fig. [2] 
for the shortest wavevectors ~ 1/L, with L the system 
size. The phase diagram is constructed again by assuming 
that the system will relax to linearly stable states charac- 
terized by the fastest relaxation rates. To linear order, the 
coupling to motor fluctuations yields oscillatory or propa- 
gating (as opposed to purely diffusive) fluctuations. Some 
of these oscillatory states are linearly unstable in a region 
of parameters, as shown in the phase diagram in Fig. [2j 

Next we consider the effect of nonlinearities by first ex- 
amining a minimal model that only retains the longest 
wavelength mode in the Fourier expansion of the nonlin- 
ear equations, Eqs. (2[3), and then comparing the latter 
to the numerical solution of the full nonlinear partial dif- 
ferential equations. 

One-mode model. — We begin by incorporating 
only the nonlinearities in the unbinding rate, while ne- 
glecting convective, elastic and pressure nonlinearities. 
We impose boundary conditions [S^mJ^^o = [d x u] x =L = 
and [(/>]x=o = [<AU=l = 0, and seek a solution of 
Eqs. ( 2j3 1 in the form of a Fourier series as, u(x,t) — 

E" =0 u m(t) cos (rhx) and <j>(x, t) = Y,m=l ^™ W sin (^)> 
where, m = mn/L. We perform a 1-mode Galerkin trun- 
cation, and only consider the dynamics of the first non- 
trivial mode, m = 1 (setting u m — <p m = 0, V m ^ 1). 
This corresponds to approximating the system by only 
its longest wavelength excitations, ignoring all the shorter 
wavelength modes. 

The coupled equations for u\ 



TTUl 



and 4>\ are given by 
B a ir 2 Ul + C2A/MT01 , (9a) 

■ (9b) 



-eui 
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These equations can be recast into an effective second 
order differential equation for u\ that has the structure 
of the equation for a Van dcr Pol oscillator |;22 coupled 
to a nonlinear spring, of the form u\ + iii [A — /(«i)] + 
Ul B a [I - f{ Ul )} = 0, with f( Ul ) = 
friction A = Tr 2 (B a - ( 2 A^) + 1. 



and 

Equations (9a) and 
([9b]) admit one fixed point (ui,<j>) — (0,0). Linear sta- 
bility analysis about this null fixed point shows that the 
fixed point is unstable (a repeller) when A < and is 
stable for positive A. From a global analysis 22 of Eqs. 



( 9a ) and ( 9b ) with the e nonlinearity, we find that the 



existence of the unstable fixed point signals the appear- 
ance a stable limit cycle as A crosses zero. In other words 
the system undergoes a supercritical Hopf bifurcation at 
A/i = A fi c i = (B + 7r~ 2 )/(£i -I- £ 2 ), with sustained oscil- 



lations for Afj, > A/i c i and stable spirals or nodes 22 
the two-dimensional, (ui,4>i) phase space, otherwise. 

We now consider the role of pressure and elastic nonlin- 
earities, while letting e = 0. Within the one-mode model 
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rived in Ref. 10 . The bifurcation diagram for the corn- 



Fig. 3: (color online) Left frame : Time evolution of ui(t) for 
various values of activity : A/i — 1.0 (black), corresponding 
to an unstrained steady state, A/i = 2.0 (red) and A/i = 3.5 
(green), corresponding to the region of sustained oscillations 
and A/i = 4.5 (blue), corresponding to the contracted steady 
state. Right frame: bifurcation diagram obtained from the 1- 
mode model. The figure shows a plot of the value uf of m at 
long times (t = 100) vs A/i. A supercritical Hopf bifurcation 
occurs at A/i = A/i c i ~ 1.5. As one increases A/i within the 
range A/i c i < A/i < A/i C 2, the amplitude of oscillations grows 
continuously. At A/i = A/i C 2 the limit cycle disappears and the 
system settles into a contracted steady state for A/i > A/t C 2. 
The dashed line indicates an unstable fixed point. Parameters: 
B — 3.0, e = 0.1. Other parameter values are same as in Fig. [2] 

we find that in this case if /3 a = a a — and B a > there 
is again only one fixed point at (ui,(j>) = (0,0) and the 
linear stability analysis is identical to that of the model 
with only rate nonlinearity. In other words, if e = 0, 
the pressure nonlinearities do not stabilize the system for 
A < 0. When a a ^ 0, (i a ^ and B a < 16a 2 /27 '/3 a 7r 2 , 
we have two new nonzero fixed points (u\,<p) = (u±,0), 
describing a contracted and an expanded state of the gel, 
respectively, with 



u± 



5a a 



± 



97r 2 /3 a 3f3 a n 



16a 2 



3tt 



2 » a B a 



(10) 



Since 5p ~ TtpQUi, u + corresponds to a contracted state 
whereas u_ corresponds to an expanded state when B > 
(A/i. The contracted/expanded state is linearly unstable 

—+(2 3 — ^2 > . 



when A!j_ = — B a 



3 4 

For A 5- < the nonzero fixed points are unstable if 
B a > 16a 2 a /27(3 a ir 2 . When B a < 16a 2 /27/3 a 7r 2 and 
Aj_ < the nonzero fixed points are stable and the or- 
bits in the two-dimensional, (ui , (f>i ) phase-space describe 
nodes or spirals settling at long times to (u±,0) 22]. 
When A± > 0, the contracted/expanded states are lin- 
early unstable. The nonlinearities in the active pres- 
sure stabilize these unstable states into stable asymmetric 
limit cycles. A supercritical Hopf bifurcation occurs when 
A/i > A/i c i = (B + tt~ 2 )/{(i + (2), which terminates to 
a stable contracted steady state for A/i > A/i c2 (A/i c2 is 
determined by the solution of \±{B, A/i) = 0). 

It is instructive to note that when a a — j3 a — 0, our 
1-mode model corresponds to that of a half-sarcomere de- 



plete 1-mode model is shown in Fig. [3j with all the nonlin- 
earities incorporated. The diagram summarizes the steady 
state crossovers of the system as we increase A/i keeping B 
fixed. Finally, the steady states of the nonlinear gel within 
the 1-mode approximation are shown in the phase dia- 
gram in Fig. [4] (indicated by dashed lines), in the (B, A/i) 
plane. The unstrained state of the gel is stable when A > 
whereas the contracted state is stable for A!J_ > 0. When 
A < and A^ < the gel exhibits sustained oscillations. 

Numerical analysis of the Continuum Nonlinear 
Model. — In this section we discuss the numerical so- 
lution of the nonlinear PDEs given in Eqs. ( 5a|5b | (but 
with no convective nonlinearities) with boundary condi- 
tions [d x u(t)] x=Q L = and [^Ja—o £ = and an initial 
strained state. We spatially discretize the PDEs using 
finite difference method and then integrate the resulting 
coupled ODEs using the rkf45 method. 

The behavior of the system as we vary B and A/i, keep- 
ing all other parameters constants, is summarized in the 
numerically constructed phase diagram, shown in Fig. [4j 
Assuming (1 = (2 = (, f ne system settles into an un- 
strained state for B > B c \ ~ 2(A/i and Hopf bifurcates 
to sustained oscillations for B < B cl . We note that the 
phase boundary in the 1-mode model B = 2(A/i — 1/7T 2 , 
although different from the numerical phase boundary 
B = -B c i, lies within the error bar of the numerics. For 
B < B C 2 ~ 0.85(A/i the system settles into a stable con- 
tracted state. The basin of attraction of the expanded 
state is much smaller than that of the contracted state. 
This is due to the choice of the sign of the coupling con- 
stants, a a and £4. The regimes predicted by our con- 
tinuum phenomenological model may be used to clas- 
sify the behaviour seen in a number of experimental sys- 
tems [5|[6j[T3]. To estimate the parameter in real sys- 
tems we assume T ~ t?/£ 2 ; with rj the viscosity of the 



permeating fluid and £ the gel mesh size 18 and use ex- 
perimental values of B. For muscle fibers, using B ~ 2 
2 x 10~ 3 Pa s |], £ = 0.5 fim, L - 100 pm 



kPa 
and k, 



23 



V 







40 ms, we estimate the dimensionless modulus 
B = B/TL 2 ko as B ~ 1. For isotropic cross-linked acto- 
myosin gels only direct measurements of the low frequency 
shear modulus G are available, with G ~ 1 — 10Pa pi. It 
has been argued, however, that these networks may sup- 
port much higher compressional forces, of the order of 
buckling forces on the scale of the mesh-size, yielding a 
value of B comparable to that of muscle fibers 24 . For 
B ~ 1 — 10 3 Pa and 77 comparable to that of water, we 
obtain B ~ 10 -3 — 1. The active coupling (A/i can be es- 
timated as [lO] (A/i ~ t;ribk m A m , where rib is the number 
density of bound motors, k m is the stiffness of the myosin 
filaments and A TO is the steady state stretch of the bound 
motor tails. Using rib ~ 10 2 — 10 4 /im -3 , k m = ApN/nm 
and A m ~ 1 ram, we obtain (A/i ~ 0.1 — 10 in dimen- 
sionless units. In the phase diagram in Fig. |4j these pa- 
rameter values would put muscle fibers in the regions U, 
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Fig. 4: (Color online) Phase diagram in the (B, A/x) plane 
for the continuum nonlinear model described by Eqs. (5a 
and (5b |. U : Unstrained, SO : Sustained Oscillations, C 
Contracted. The points are obtained by numerical solution 
of the full PDE's; the error bars are determined by the step 
size used in the A/i increments. The dashed lines indicate the 
1-mode phase boundaries. Inset : Plot of extension A£(t) — 
u(0,t) - u(L,t) for B = 3.0 and Afi = 1.0 (black), Afi = 2.0 
(red), A/i = 3.5 (green) and Afi — 4.5 (blue). Parameter 
values are same as in Fig. [3] 



SO and C respectively as we increase n b . Isotropic ac- 
tomyosin networks may have a much lower value of B, 
possibly corresponding to the contracted region. Indeed, 
while spontaneous contractility has been observed in vitro 
in reconstituted actomyosin networks [5j , no experimental 
evidence has yet been put forward of spontaneous oscilla- 
tions in these systems. 

Conclusion. — We have presented a generic contin- 
uum model of a cross-linked active gel which can be used 
to describe a wide variety of isotropic elastic active sys- 
tems. We find an elastic active gel system can, in gen- 
eral, be tuned through three classes of dynamical states 
by increasing motor activity: an unstrained steady state 
of homogeneous constant density, a state where the lo- 
cal density exhibits sustained oscillations, and a sponta- 
neously contracted steady state, with a uniform mean den- 
sity. The continuum model is not strictly hydrodynamic 
due to presence of the fast variable <j), describing the dy- 
namics of motor proteins. The motor binding/unbinding 
kinetics plays a crucial role in generating oscillating states 
at finite wavevectors. The one-mode model is in excel- 
lent qualitative agreement with the results described by 
the solutions to the nonlinear continuum equations and is 
comparable to a variety of one dimensional models for ac- 
tive elastic systems. Quantitative agreement in the phase 
boundaries between the one-mode model and the contin- 
uum model fails due to the non hydrodynamic nature of 



the model and we expect the oscillatory and/or contractile 
behaviour to depend on system size [25]. 
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